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We report on tensor reduction of five point integrals needed for the 
evaluation of loop-by-loop corrections to Bhabha scattering. As an exam- 
ple we demonstrate the calculation of the rank two tensor integral with 
cancellation of the spurious Gram determinant in the denominator. The 
reduction scheme is worked out for arbitrary five point processes. 

PACS numbers: 13.20.Ds, 13.66.De 



1. Introduction 

Evaluation of the Bhabha cross-section with two-loop accuracy [1~4] 
would be incomplete without taking into account contributions from one 
loop diagrams with photon emission from internal electron lines. So far 
bremsstrahlung from external legs (lowest order) has been taken into ac- 
count [5] and implemented in [6-8]. 

The general form of the five point tensor integrals in the diagrams of 
Fig. 00 is: 

J in'^/^ C1C2C3C4C5 

where Ci = {q + qiY — mf, i = 1, . . . , 5. We also make the (arbitrary) 
choice ^5 = 0. 
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Diagram 1 01 ; topology u95_ 
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Diagram 102. topology u96_ 



Diagram 103; topology U97_ 





agram 104; topology u98_ 



Diagram 1 05; topology u99_ 




Diagram 106; topology uiO0_ 

— ' — wvvwvWaa/w* 



n 107; topology u101_ 




agram 108; topology u102_ 




Fig. 1. Eight diagrams at one loop level obtained using DIANA [11]. 



Massive five point functions tend to be unstable numerically in certain kine- 
matic domains due to the appearance of an inverse Gram determinant. Cal- 
culating the tensor integrals directly, e.g. by using Mellin-Barnes (MB-) rep- 
resentations, we have to cope with five dimensional MB-integrals (before ex- 
pansion in epsilon). After e expansion, five dimensional MB representation 
still remains at (e.g. for second order tensor). The MB-reprcscntations 
are derived with the aid of the Mathematica packages AMBRE and MB [9, 10]. 
They are hard to evaluate both numerically (and in fact presently the MB 
package works in Euclidean region only) and analytically. A reduction of 
tensor five point functions to scalar four and three point functions solves 
that problem. 



2. Reduction of five point functions 



We give examples of reduction for the scalar and second rank tensor 
five point function where the spurious Gram determinant in denominators 

is avoided (sec also [12]). 

To make a reduction one can follow the standard Passarino-Veltman [13] 
reduction scheme, or, as we will do, follow a scheme based on [14-16]. For a 
diagram with internal lines 1 ... n we can introduce the so called "Modified 
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Cayley Determinant": 
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(2.1) 



with Yij = —{qi — Qj)^ + in^ + m^. By cutting from ()5 rows ji, j2, • • • and 



columns ki,k2, ■ ■ ■ we get so-called " Minors" . The sign of a " Signed Minor" , 

(2.2) 



Ji 32 ■■ 
ki k2 . . 

is determined by the sum of indices of excluded rows and columns and by 
taking into account the appropriate signatures of the permutations, taken 
separately from excluded rows and columns (important when cancelling 
more than one row and one column), e.g. signatures of permutations for 
(54) is: + for rows and — for columns. In this way: 



'0^ 
,0, 



Yn Y12 
Y12 Y22 



Ym Y2 



2n 



Ym 

Y2n 
Y 



(2.3) 



Using " signed minors" one can write tensor reduction formulas in a compact 
and elegant way. 

We start from the simplest case of a scalar five point function. As 
starting point we use the recursion relation [15]: 



(d-^., + l)G„_,4'^+^)= 2A. + ^(|^) 

i=l L k=l \ ^ / 



lid) 
^n ' 



where 



A. 



'0^ 
,0, 



Gn- 



E(l^l=2()n. 

k=l 



, ft. . 



(2.4) 



(2.5) 



The operator k~ decreases the A;-th propagator power by one. Because in 
our case powers of propagators are all equal to one, and we are interested 
in the case with n = 5, the above relation changes into: 



(^-4)()5/^' 



^0^ 
.0. 



5 

E 

k=l 



rk 
J4 . 



(2.6) 
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In the limit d ^ 4 we get as final reduction formula for the scalar five point 
function: 

^5 = 7^E(°)/|, (2.7) 
W5 k=i VVs 

where is the scalar four point integral obtained by cancelling the A;-th 
line in the five point function diagram. 

For the vector integral we refer to [15]. The second rank tensor integral can 
be written in the form: 

ir=j:ifl^n.A'f-l9^^ltK (2.8) 

Here [d+Y = 4 + 2Z — 2e, and = 1 + Sij. We start from the following 
formulas [15]: 

„ . . T[d+? _ ©5 Tld+] , V ^^^5 Ad+hs , Ad+] q\ 

W5 s=i U5 U5 



jj=l 1/5 



,-[<i+],s _ \is)5 Ts I W5 Tst to ^^\ 



Using all the above contributions we end up with: 



I^"" = E (2-12) 

2J = 1 



where 



T _ J \^ r^^M r« VjVsVis^S I Vj/ 5 Vis/ 5 jst 



(2.13) 

However, we can avoid the determinant Os in the denominator. We begin 
with the following structure: 

if = [/r - Eoogn + Eoog^'' (2.14) 
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The problem of finding an appropriate ansatz for £"00 fias been solved in [16]: 

5 



00 



\0) s=l 



00' 



(2.15) 



where Dqq is the g^'^ term of the tensor four point function. The complete 
result for Eqq is: 



1 1 



00 



t=i 



, ts . 



jst 
^3 



(2.16) 



To demonstrate in details how Os is cancelled in the square bracket in (j2.14p . 
we have to consider four and three point functions separately by analyzing 
the coefficients of I| and /J* in (j2.9p . subtracting the g^" term according to 
(I2T0]) and (|2T6]) . For /| we have 



1 



si 5 



()505( 



vs 



, is , 



+ 



,0s. 
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(2.17) 



i.e. we have to show that indeed Os cancels and we have to give an explicit 
expression for Xfy 

A useful property of Xfj is its symmetry w.r.t. the indices i and j for 
fixed s. Obviously the third term in the curly bracket of (j2.17p is symmetric 
since we consider a symmetric determinant. The symmetry of the first two 
terms means 



+ 



s ' 
,0, 



5 I 
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(2.18) 



The first square bracket of (|2.18p can be evaluated using (A. 13) of [14], i.e. 

(2.19) 



'0^ 
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and (I2.18P then results in: 

This is proved by multipUcation with Os and using (A. 8) of [14] with r = 

Inserting this, products of three factors of the form (^)g cancel pairwise, 
q.e.d. . 

Further, the following relations (A. 11) and (A. 12) of [14] are important. 

As simplest case we now immediately obtain from (I2.17P {• • - j^^ = 0, i.e. 
Xl^ = 0. Applying (12:221) to (12171) . we see: 

and due to the symmetry in i and j we also have: 

which gives us a hint of how Xfj might look, namely due to (j2.23p it should 
contain a term A further contribution, summed over, must 

vanish. Due to Xf„ = it must contain a factor (^■^■)^. The second factor of 
this further contribution can only depend on s and has been determined by 

Os/5' 



explicit calculation to be ^ ■ Thus we conclude 



For /I* a slight generalisation yields the coefficient 
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with X^: = Xff ^ and the final result is: 

€= E Q^<l3h,^J+9'"'Eoo, (2.27) 

where 

1^1 ^ 

■^5,ij = Tor E - J2 ^ij^s) ■ (2-28) 

lo/S s=l f=i 

Xj?* is defined in (j2.26p and Eqq in l2.16i /J* is the scalar three point function 
obtained by cutting the s and t lines in the five point diagram. 

3. Automatization 

We have developed a Mathematica package which reduces general tensor 
five point functions up to rank three. Using LoopTools [17] we have per- 
formed numerical cross-checks which ensured us about the correctness of the 
package. We have also made independent checks with sector decomposition 
and the Mellin-Barnes method. Some examples are given in [18]. 

In summary, we have developed tools to deal with five point Bhabha 
amplitudes. The next step is to implement it into the complete numerical 
calculation. 
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